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ABSIRACT 


freseoseCeeOr TETS “STUEY™ WAS~ TG DEVELOP: A HIGHLY 
ACCURATE EGNG RANGE CECMPUTER PROGRAM TO. DOO PREDICTIONS OF 
TROPOSPHERIC RADIO TRANSMISSION LOSS OVER IRREGULAR TERRAIN. 
fae COMPUTER PROGRAM WAS TO BE EASY TO USE WITH AN INPUT: AND 
OUTPUT DATA FORMAT SUCH = THAT THE USER CELLD-USE TRE PRGGRAM 
Wetm eh ttTEe PREPOR TNSTRUCTICN. SIGNAL STRENGTH CALCULATION 
Sepeeteiiy WAS ALSO TC BE INCLUCEC IN TRE PROGRAM. THE BASIS 
Bere inie PROGRAM TREORY WAS ESSA RESEARCE REPORT NUMBER ERL 


Seats 6f, “PREETCT ICN OF TROPESPHERIC RADIO TRANSMISSION 


LOSS OVER IRREGULAR TERRAIN." 
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I. INTROCLCTION 


WHEN RF POWER RACIATES FRCM A TRANSMITTING ANTENNA IT 
IS SPREAD OVER A WIDE AREA. THUS, THE PCGWER AVAILABLE AT THE 
RECEIVING ANTENNAS IS CNLY A SMALL PART OF THE RADIATED 
POWER. 

THE LOSS DUE TO TRANSMISSICN BETWEEN THE TRANSMITTING 
AND RECEIVING ANTENNAS DETERMINES IF THE RECEIVED SIGNAL 
WILL BE USEFUL. PCOR GUALITY CR PCOR RELIABILITY CAN RESULT 
IF THE MAXIMUM ALLCWABLE TRANSMISSION LOSS FOR THE RADIO 
SYSTEM IS EXCEEDED. CN PATHS THAT APPROXIMATE THE IDEAL FREE 
SPACE OR PLANE EARTH, FAIRLY GOOD PREDICTIONS OF 
TRANSMISSION LOSS CAN BE MADE. THE ACTUAL PATH GECMETRY OR 
ATMOSPHERIC CONDITIONS IN THE REAL WORLD CIEFER SC MUCK FROM 
THE BASIC ASSUMPTIONS THAT ABSCLUTE ACCURACY CANNOT BE 
EXPECTED. 

THERE IS GREAT INTEREST TODAY IN. RACIO WAVE 
PROPAGATION BEHAVIOR IN THE TROPCSPHERE; WHICH IS THE LAYER 
OF THE ATMOSPHERE WHERE CLOUDS CAN EXIST. ITS UPPER LIMIT; 
MMEEDN THE TROPOPAUSE, IS AT A HEIGHT OF ABCLT SIX 
KILOMETERS AT THE POLES AND ABCLT EIGHTEEN KILCMETERS AT THE 


EQUATOR. 


MANY COMPLEX FCRMULAS HAVE BEEN CERIVED TO MAKE BETTER 
PREDICTICNS OF TROPOSPKERIC RACIO TRANSMISSION LOSS OVER 
IRREGULAR TERRAIN. THESE EQUATIONS HAVE BEEN PUT TOGETHER IN 
ORDER TO ESTABLISH A PROPAGATICN THEORY. THE THEORY HAS BEEN 
TESTED AGAINST A LARGE NUMBER CF PROPAGATICN MEASUREMENTS 
AND HAS HELD UP QUITE WELL. THE MORE ACCURATE THE THEORY THE 
MORE EQUATIONS NEEDED. 

IN THIS DAY ANC AGE ACCURATE THEORY REQUIRES A LOT OF 

\ 
TIME TO SOLVE MANY CCMPLEX EQLATIONS. A COMPUTER IS NEEDED 
TO COMBINE THE NUMERGUS EQUATICNS AND SOLVE THEM IN A SHORT 
AMOUNT OF TIME GIVEN CERTAIN NEEDED INFGRMATICN. BEFCRE A 
COMPUTER CAN BE USED IT MUST HAVE A SOFTKARE PRCCRAM WHICH 
COMBINES THESE WELL ESTABLISHEC EQUATICNS IN AN ORDERLY 
FASHION IN ORDER TO ARRIVE AT A NUMERICAL SOLUTION. 

IN JULY 1968 ESSA RESEARCE LABORATORIES PUBLISHED A 
TECHNICAL REPORT TITLED: PRELICTION OF TROPOSPHERIC RADIO 
: TRANSMISSICN LOSS OVER IRREGULAR _TERRATN-_A_ COMPUTER METHODS 
1968 BY A. Ge LONGLEY AND P. Les RICE. THE REPORT NUMBER WAS 
ERE 79-I1TS 67. 

THE ESSA REPORT CONTAINED A COMPUTER METHOD FOR 
PREDICTING LONG TERM MEDIAN TRANSMISSION LOSS OVER IRREGULAR 
TERRAIN. THE METHOD WAS APPLICABLE FGR RACIG FREQUENCIES 
ABOVE 20 MHZ AND MAY BE USED EITHER WITH CETAILED TERRAIN 


PeGenees, IF ORee ACTUAL “PATHS CR WITH PROG DEES “THAT ARE 
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REPRESENTATIVE OF MECIAN TERRAIN CKARACTERISTICS FOR A GIVEN 
AREA. ESTIMATES OF VARIABILITY IN TIME AND WITH LOCATION AND 
A METHOD FOR COMPUTING SERVICE PROBABILITY WERE INCLUDED. 

Tie CO ROGnAM.: GiVEN = IIN’ THE REPORT -WAS° VERY LIMITED IN 
Beas TRECT APPLICATIGNS.s IT WAS WRITTEN FOR USE CN A CCNTROL 
DATA CORPORATICN COMPUTER, WHICrF MADE IT UNUSEABLE FOR AN 
PoMeCOMPUTER. THE PRGCRAM, AS WRITTEN, COULD NOT READ IN 
Sera PREM A CARD SOURCE. THE WORK INVOLVED IN THIS THESIS 
Osts NO NEW PROPAGATION THEORY, HOWEVER, [T DOES COMPLETELY 
REWORK THE SOFTWARE COMPUTER PROGRAM IN THE REPCRT IN ORDER 
BeeMAKE IT MORE EASY 710 USE, MCRE FLEXIBLE ANC ACCURATE, AND 
PeeUDES IN THE GUTPUT, A GRAPHICAL PLCT. IT ALSC RETURNS 
SIGNAL STRENGTH. 

THE PROPAGATION PROGRAM IS INTENDED FOR USE WITHIN THE 


FOLLOWING RANGES: 


PARAMETER RANCE 
FREQUENCY 20 TO 40,000 MHZ 
ANTENNA HEIGHTS Oso 18 34000 M 
DISTANCE 1 TO 2,000 KM 
SURFACE REFRACTIVITY 250 T0 400 N-UNITS 
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hrmesANGLE OF ELEVATION OF EACH HORIZON RAY ABOVE THE 
HORIZONTAL SHOULD NOT EXCEED 12 DEGREES, AND THE DISTANCE 
FROM EACH ANTENNA TO ITS HORIZCN SHOULD NOT BE LESS THAN ONE 
DeWtis OR MORE THAN THREE TIMES, THE CORRESPONDING 


SMOOTH-EARTH DISTANCE. 
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II. COMPUTER FROGRAMMEC TRANSMISSION THEORY 


BEFORE GGING INTC THE PARTICULAR ASPECTS GF HOW TO USE 
THE PROGRAM, A BRIEF GENERAL CUTLINE OF THE THEORY USED IN 
THE ERL 79 — ITS 67 REPORT WILL BE GIVEN. 

RADIO TRANSMISSICN LOSS IN TROPOSPHERIC PROPAGATION 
DEPENDS ON CHARACTERISTICS CF THE ATMOSPHERE ANC TERRAIN. 
THE REFRACTIVE INDEX GRADIENT AEAR THE EARTH'S SURFACE IS 
THE MCST IMPCRTANT ATMOSPHERIC FARAMETER FOR PRECICTING A 
LONG TERM MEDIAN REFERENCE VALLE OF TRANSMISSION LOSS. THIS 
SURFACE GRADIENT LARGELY DETERMINES THE BENDING OF A RADIO 
RAY AS IT PASSES THROUGH THE ATMOSPHERE. RAYS MAY BE 
REPRESENTED AS STRAIGHT LINES, bhITHIN THE FIRST KILCMETER 
ABOVE THE EARTH'S SURFACE, IF AN “EFFECTIVE EARTH RADIUS", 
1S DEFINED AS A REFRACTIVITY. 

FOR THE CALCULATICN OF TEE LONG-TERM REFERENCE VALUE, 
THE MINIMUM MEAN VALUE 1S PICKED TO CFARACTERIZE AVERAGE 
ATMOSPHERIC CONDITIONS. THE MINIMUM MONTHLY MEAN VALUE OF 
SURFACE REFRACTIVITY I'S OBTAINED FROM MEASUREMENTS FROM SUCH 
MAPS AS BEAN, HORN, AND CGZANICH (1960), "LL CLIMATIC CHARTS 
AND DATA GF THE RACIC REFRACTIVE INCEX, FOR THE UNITED 


STATES AND THE WORLD," NBS MCACGRAPH NG. 22 (U.S. COVERNMENT 
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PRINTING OFFICE, WASHINGTON, CO. Ce) 

Sr hodenewt Ce wEERECTS,;, SUCH AS CHANGES .IN THE REFRACTIVE 
INDEX; CHANGES IN THE AMCUNT OF TURBULENCE OR 
STRATIFICATION, ABSORPTION BY CXYGEN, WATER VAPGR; CLOUCS, 
AND PRECIPITATION, ARE ALLOWED FCR BY EMPIRICAL ACJUSTMENTS. 

PROGRAM CALCULATIONS ARE TNCLUDED TC HANCLE EFFECTS OF 
MeerOsSPHERTC SCATTER AND DIFFRACTICN OUE TQ KNIFE EDGES, 
ROUGH TERRAIN, ANC SMCOTH SPHERICAL EARTH. SHORT AND LONG 
TERM FADING ADJUSTMENTS TO MEDIAN ATTENUATIONS ARE ACCOUNTED 
FOR BY DERIVING FACTORS WHICH COME FROM AVERAGING FADING 
MEASUREMENTS MADE CN MANY OLE FERENT TREPCSPHERIC 
COMMUNICATIONS SYSTEMS. 

fEeverAl CED SPREEPEES ARE AVAILABLE FOR SPECIFIC PATHS, 
TRANSMISSICN LOSS CAN BE CALCULATED. TO CHARACTERIZE 
TERRAIN, PROFILES MAY BE REAC AT REGULAR INTERVALS IN BOTH 
N-S AND E-W OIRECTIONS OR ACTUAL CR PROPOSEC TERRAIN 
CONDITIONS MAY BE COMBINED Tee eROVODE GAs eSENGLE “SET OF 
PROFILES EGR WHICH “AN ESTEMATE OF WEDIAN FRCPAGATION 
CONDITIONS IS.DESIRED. THE INTERDECILE RANGE CF TERRAIN 
felch S ABOVE AND BELCW A FIXEC LINE FITTEC TO ELEVATIONS 
PEEMERSEA LEVEL, IS CALCULATEE AT FIXED CISTANCES. USUALLY 
MEDIAN VALUES OF THE INTERDECILE RANGE INCREASE WITH PATH 
BPENGTH TO AN ASYMPTOTIC VALUE WHICH IS USED 10 CHARACTERIZE 


THE TERRAIN. 
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ieee OF TERRAIN ASTMeTOotIGC VALUE (METERS) 
SE INTERGECTEC HELGHT 


VERY SMOOTH PLAINS OR WATER 0=5 
SMCCTH PLAINS S= 2.0 
SeIGHELY ROEEING PLAINS 2C-40 
ROLLING PLAINS 4C-&0 
RES 807150 
MCUNTAINS 150-300 
RUGGED MOUNTAINS 300-700 
EXTREMELY RUGGED MCUNTAINS > 100 


POVEAECUEATE REFERENCE ‘VALUES OF TRANSMISSTCN. LESS FOR 
ReeOPEGIEIC “APPLICATION, A MINIMUM OF FCUR ESSENTIAL 
PARAMETERS MUST BE SUPPLIED. THESE  FCUR PARAMETERS ARE 
FREQUENCY IN MEGAHERTZ; * THE’ TRANSMITTING AND RECEIVING 
SN beENNA ~ *FEfCHiS ~ABECVE GROUNCE “IN METERS,  ANO’ TRE PATH 
DESTANGE IN KILOMETERS. PATH PARAMETERS SUCH AS HORIZON 
DYSWANCES AND ELEVATIGN ANGLES, MAY ALSO BE DERIVED FROM THE 
PRECEDING PARAMETERS. 


the CONDUCTIVITY OF THE EARTH®S’ SURFACE AND ITS 
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peer ivVitY OR RELATIVE DIELECTRIC CONSTANT’ ENTER INTO THE 
EeeGULATIGNS FOR  LINE-GF-SIGHT AND DIFFRACTION ATTENUATION. 
IF THE VALUES OF THESE GROUND CCNSTANTS ARE NOT KNOWN, THE 


FOLLOWING TABLE VALUES MAY BE ASSUMED: 


reer OF SURFACE CONDUCTIVITY (MHO/M) DEEL ECTRIC, CONSTANT 


POOR GROUND 0.COl 4 
AVERAGE GROUND 0.005 15 
GOOD GROUND 0.020 BS) 
SEA WATER 5-000 81 
BRESH WATER 0.010 81 


Pi esUGERCTENTIOY VEGW PREQUENCTES, THE EFFECT OF THE 
CONDUCTIVITY DOMINATES WHILE AT SUPE LEBEN TEY HIGH 
FREQUENCIES, THE DIEEECERTE CONSTANT DOMINATES. THE 
TRANSITION OCCURS BETWEEN 3200 16 3000, MEGAHERTZ FOR OVERSEA 
TRANSMISSION, WHILE CVER “AVERAGE GROUND THE TRANSITION 
OCCURS BETWEEN 5 TO 50 MEGAFERTZ. FOR - PROPAGATICN OVER 
IRREGULAR TERRAIN, AT FREQUENCIES ABOVE 100 MEGAHERTZ: WITH 
ANTENNAS MORE THAN FIVE METERS ABOVE GROUND, THE EFFECTS OF 
THE GROUND CONSTANTS ARE SLIGHT AND THE RESULT, FCR VERTICAL 


AND HORIZONTAL POLARIZATION IS NEARLY THE SAME. UNDER THESE 
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SenDITIGNS THE METECC IS SIMPLIFIED BY ASSUMING THE 
Peo UE OP THE THEORETICAL REFLECTION CGEFFICIENT EQUAL TO 
0.95 AND’ THE PHASE “SHIFT  ECUAL nee cERO.  ReR MANY 
Peeve lOvs this RESUETS “IN AN ESTIMATE OF THE EFFECTIVE 
PerPeeCtLON COEFFICIENT GF O.9, ANC THE ATTENUATICN MAY BE 
PePeccateD DIRECTLY. THESE APPREXIMATIGNS ARE NET APPLICABLE 
RoR TRANSMISSION GVER THE SEA. 

TO] SUMMARIZE, THRE PARAMETERS REQCUIREC FOR COMPUT ING 
REFERENCE VALUES OF TRANSMISSION, OR THE CORRESPONDING 
ATTENUATION BELOW FREE SPACE ARE: FREQUENCY IN MEGAHERTZ, 
DISTANCE IN KILOMETERS, ANTENNA HEIGHTS ABCVE GRCUND IN 
Meeks, ESTIMATES CF SURFACE REFRACTIVITYs THE. TERRAIN 
PeReEGUEARTTY EN METERS, THE CENDLCTIVITY IN MHOS/METER, AND 


Sao keCATT VE DIELECTRIC CONSTANT CF TRE GROUNC. 
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Ill. USING THE MODIFIED ERL PROGRAM 


Mae MODIFIED ERE COMPUTER PROGRAM WAS WRITTEN FOR EASE 
Seeeuteutr READING, ACCURACY, AND EASE OF USER DATA INPUT. 
THE PROGRAM LISTING CONTAINS INSTRUCTIONS FOR DATA’ INPUT. 
THIS ALLOWS THE PROGRAM TO BE USED WITHOUT EXTERNAL PROGRAM 
INSTRUCTIUNS.~ 

BHeEReE ARE GNEY FCUR DATA INPUT CARDS USEC IN THIS 
PROGRAM. IN MANY CASES ONLY TWO DATA INPUT CARDS ARE USED. 
fms FPF ROGRAM PREDICTS LONG-TERM MEDIAN RADIO TRANSMISSION 
BeooeOVER TRREGULAR TERRAIN. THUS, THE USER  SKCULD NOT 
Seer A SHORT TERM “LOSS PREDICTION GQUTPUT EVEN THOUGH THE 
PNPUT DATA IS OF A -SHCRT TERM NATURE. 

INTEGER NUMBERS MUST BE FUNCHED SO THAT THE NUMBER 
ENDS IN THE LAST COLUMN NUMBER FIELD. TE SPACES ARE EERT JAm 
THE END OF THE FIELD THEY WILL BE READ AS ZEROS; PENS ERECT 
MULTIPLYING THE DESIRED NUMBER BY A POWER OF TENe REAL 
NUMBERS ARE PUNCHED AS) 7A STRING “OF DIGITS CONTAINING A 


DECIMAL, AND MAY BE PUNCHED ANYWEERE IN THEIR FIELD. 
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FIRST DATA CARD: 
Pies seikols DATA CARD. CENTAINS CATA’ -COVERING CODE 
SuMoens FOR TYPE OF TERRAIN, CE LOSS INFORMATION DESIRED, 


AND DISTANCE BETWEEN ANTENNAS. TEFIS CARD MUST BE INCLUDED IN 
tie DATA CARD SET. 


CobUNN=: |Ol—-10 — “TYRE OF TERRAIN’ -INTEGER 


CODES I COLCRADO PLAINS 
COMTT-TERIRD: DATA CARD) 
CODE: -2 COLOKRACO MOUNTAINS 
; (OMIT THIRD DATA CARD) 
CODE: 3 OEE 
(OMIT THIRD DATA CARD) 
CODE: 4 USER TERRAIN CATA 
WITHOUT DATA SUPPRESSION 
CODE: 5 USER TERRAIN CATA 


NO “OUTPUT PARAMETERS” PRINTED 
COLUMN=s Pi—2o "De LOSS" -—INTEGER 


CEDE= 0 NG Se IATA DES TRED 
(OMIT SECOND CATA CARD) 
CODE 1 DB DATA DESTREC 


COLUMN: 21-30 “CISTANCE BETWEEN ANTENNAS" ~REAL 
(1-2,0C0) K~METERS 
SECOND DATA CARD: 

THE SECOND DATA~- CARD CCNTAINS DATA COVERING SIGNAL 
STRENGTH DB CALCLLATICNS AND CLTFUT. IF COLUMN 11-20 CF THE 
FIRST DATA CARD CONTAINS A “O" THE SECOND DATA CARD IS 
OMITTED. | 

COLUMN: O1-10 "TX POWER CLT" (WATTS) 

POSITIVE REAL NUMBER 

COLUMN: 11-20 “TX ANTENNA GAIN" (C2) 

POSITIVE REAL NUMBER 
~ COLUMN: 21-30 "RX ANTENNA GAIN" (08) 
POSITIVE REAL NUMBER 
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COLUMN: 31-40 "TRANSMITTER PINE LOSS. Cb B) 
POSITIVE REAL NUMBER 
CEPUMN: 41750 “REQUIR'IEE RECEIVER LINE LOSS" (DB) 
POSITIVE REAL NUMBER 
COLUMN] 51-60. “RECEIVER INPUT LEVEL™ (CBM) 
POSITIVE REAL NUMBER 
THIRD DATA CARD: 

SHE TALRD DATA ‘CARD CONTAINS | DATA ‘COVERING SURFACE 
BeERACTIVITY: SURFACE CONDUCTIVITY, PERMITTIVITY CR RELATIVE 
PPeeectTRIC CONSTANT, AND INTERCECILE RANGE. TEIS CARD IS NOT 
Seo el COLUMN O1-10 CF THE FIRST DATA CARD CONTAINS A "1", 
ees OR M3, 

CORUMN=s OL-=10 “SURFACE REFRACTIVITY" -INTECER 

(250=4C0) (N= UNITS -NS 
SGEUMN=  Pi-~20 "SURFACE CONDUCTIVITY” (MHGS/ METER) 
COECUMN= 20-305 “PERMITTIVITY GR RELATIVE DIELECTRIC 
CON STAN GO 
CGEUMN=s 31=40 “TNTERDECILE RANGE" 
DEETASK 
FORTH DATA CARD: 

THE FORTH DATA CARD CONTAINS DATA COVERING ANTENNA 
POLARIZATION, FREQUENCY, AND ANTENNA HETChT OF THE 
TRANSMITTING AND RECEIVING ANTENNA’ THIS CAREC MUST BE 
BNCEUDED IN THE CATA CARD SET. 

CCLUMN: 01-10 “ANTENNA FCLARIZATICN® +~REAL 

CODE: 01.00 © VERTICAL POLARIZATION 
CODE = 1500 HORIZCNTAL PCLARIZATICN 
COLUMN= LI—-20 “FREQUENCY (MEGAKERT Z) 
(20-4C,CCC) 
COLUMN: 21-30 "STRUCTURAL ANTENNA FEICGHT CF 
TRANSMITTING ANTENNA" ~—REAL 
(Q55-2C00.0) METERS 
COLUMN:31-40 "STRUCTURAL ANTENNA FEICHT OF 
RECEIVING ANTENNA® -REAL 


COl.5-2C00.20) METERS 
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DATA CARD GONE “CONTAINS —CCCES FOR TRE USE OF COLORADO 
PLAINS; COLCRADQ MOUNTAINS, GR OHIO TERRAIN PARAMETERS. 
MHESE PARAMETER CPT IONS WERE sPUT “IN THE PROGRAM SINCE’ THEY 
Seen WELE-USED GENERAL TYPES CF DIFFERING TERRAIN. FOR 
Beeamele, IF THE USER DOES NOT HAVE ANY PARTICULAR — TERRAIN 
CONSTANTS TO USE AND WANTED TC CHECK HIS DATA USING A 
MOUNTAINOUS TERRAIN, HE WOULD USE THE COCE IN TEE FIRST DATA 
CARD FOR COLORADQ MCQUNTAIN TERRAIN. THE FOLLOWING TABLE 
fev es ThE CONSTANTS USED -IN TEE PROGRAM FOR THE PRESET 


TERRAINS: 


SURFACE SURFACE DIELECTRIC DH 
WERRADN REFRACTIVITY CONCUCTIVITY ~—CONSTANT CELTA—-H 


OHIO Sr250 0.CC5 15.0 90.0 
COLGRADC PL. 290.0 0.CC5 Lo 20 90.0 
COLORADO MT. 290.0 0.0C5 15.6 €50.0 


EXAMPLE PROGRAM DATA CARD SETUP 
THE FOLLOWING TG. “EXAMPEES SHGULE GIVE ) ThE USER A 
BemieR GRASP GF SETTING UP TEE CATA CARDS FOR . TRE MODIFIED 


ERL PROGRAM. 


EXAMPLE ONE: 
THE . USER WANTS -TO USE..CCLCRADO PLAINS. TYFE TERRAIN 
DatAe HE ALSO WANTS TG USE THE FOLLOWING DATA: 
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ee ee TES 
Beg ueENC Y 9 =—=——— rn ne oe ne 450.200 MEG-HERTZ 
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NOTE: 


AEE INTEGER 


NeES EXAMPLE (CATA CARED. SETTINGS 
EXAMPLE 
GNE 


— eye ee ae ee ee 


1.C 

45C.C 

4.C 

4C.C 

VADEES MUST BE 


RESPECTIVE CELUMN FIELECS.« 


TRE NEX® SECTICN,s 
SAFPEE CUTPLT DATA FCRMATS. 


PAGE 23. 
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